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Abstract 

We apply our technique on the dispersion of magnetic fields in molecular clouds to 
high spatial resolution Submillimeter Array polarization data obtained for Orion KL in 
OMC-1, IRAS 16293, and NGC 1333 IRAS 4A. We show how one can take advantage of 
such high resolution data to characterize the magnetized turbulence power spectrum in 
the inertial and dissipation ranges. For Orion KL we determine that in the inertial range 
the spectrum can be approximately fitted with a power law fc-(2-9±0.9) 

we report a 

value of 9.9 mpc for AaDi the high spatial frequency cutoff presumably due to turbulent 
ambipolar diffusion. For the same parameters we have ~ /j.-(i-4±0.4) ^^^^ tentative 
value of Aad ^ 2.2 mpc for NGC 1333 IRAS 4A, and ~ k-^i-S±o.3) ^j^j^ ^^^^j. ^j^^j^ 
of Aad ^ 1-8 mpc for IRAS 16293. We also discuss the application of the technique to 
interferometry measurements and the effects of the inherent spatial filtering process on 
the interpretation of the results. 

Subject headings: ISM: clouds — ISM: magnetic fields — polarization — turbulence 



1. Introduction 



The fact that the relative import ance of magnetic fields and t urbulence in the process of star 
form ation is still a matter of debate (jMouschovias &: TassisI l2009l : ICrutcher. Hakobian. &: Troland 
20101 ) can be traced to the many difficulties of probing magnetic fields in molecular clouds. While 
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the Zeeman effect still provides the only direct way of measuring the strength of (generally the line 



of sight component o f) magnetic fields (jHeileslll997l : ICrutcher et al.lll999l : iBrogan &: Trolandll2001 



Falgarone et al.ll2008l ). its weakness in the interstellar medium limits the types of environments and 
number of regions where detections can successfully be obtained. Without the existence of a prolific 
technique for measuring the strength of all components of the magnetic field vector, it is difficult to 
precisely quantify its importance. It is therefore imperative to keep seeking new ideas and techniques 
that could provide information, even partially, on the nature of magnetic fields. Interestingly, such 
newly introduced methods of observation or analyses that have been recently put forth for the study 



Li & Houde 


2008; 


Hever et al. 


2008) 



In this paper we continue our previous w ork on the characterizati on of magnetized turbulence in 
molecular clouds using polarization maps. In Hildebrand et al. (j2009l ) (hereafter Paper I) we showed 
how the turbulent to ordered magnetic field strength ratio can be evaluated through the structure 
function of the polarization angles (i.e., the dispersion func tion) without assum ing any model for the 
ordered component of the magnetic field. Subsequently in llloude et al.l ( 2009 ) (hereafter Paper II) 
we generalized this analysis by including the process of signal integration through the thickness of the 
cloud and across the telescope beam. An important development of Paper II was the determination 
of the magnetized turbulent correlation length scale, which is in effect a measure for the width of 
the magnetized turbulence power spectrum. We go further in this paper by taking advantage of 
high resolution polarization maps to characterize the magnetized turbulence power spectrum in the 
inertial and dissipation ranges. 

We will start in Section [2] with a brief review of the material presented in Paper II on the 
cloud- and beam-integrated dispersion function that is necessary for our analysis. In Section |3] we 
present the data on which we will perform our analysis: high resolution SMA polarizatio n data of 
Orion KL (similar to those used for the map presented in Fi g. 2 c) of iTang et al.ll2010l ). as well 
as the previously pubhshed SMA data of NGC 1333 IRAS 4A (joirart et ahlbood ) and IRAS 16293 
(jRao et al.ll2009l ). We present the results of our dispersion function analysis in Section HI where we 
emphasize some important considerations that are specific to interferometry data. We then follow 
up with a discussion on the limitations of our technique, improvements that could be implemented, 
and future applications in Section [H We end with a short summary in Section [H Finally, an 
appendix pertaining to some data processing aspects can be found at the end of the paper. 



2. The Magnetized Turbulence Power Spectrum 

Following the analysis of Paper I and Paper II we model the dispersion of the difference A$ (i) = 
$ (x) — <I> (x + £) in the polarization angle <I> measured at two positions separated by a distance £ 
with 
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(cos [A^ {£)]) 



(B-B(0)> 



where (• • • ) denotes an average and £ = \£\ and (^B-B(£)y 
cloud- and beam-integrated magnetic field is defined with 



(B (r) •B(r + £)) (see below). The 



B(r) 



1 

- F (a, z) B (a, z) dz 



d a, 



(2) 



where the beam profile is denoted by H (r), while the weighting function F (a, z) > is the polarized 
emission associated with the magnetic field B (a, z) and A is the maximum depth of the cloud along 
any line of sight. The quantity r is the two-dimensional polar radius vector on the plane-of-the-sky 
and z the depth within the cloud. That is, the position vector in the cloud is given by 



X = rer + zBz (3) 

with e.f. and the unit basis vectors along r and the z-axis (which is oriented along the line of sight), 
respectively. The distance i in Equation ([TJ is also confined to the plane of the sky. We assume 
that the magnetic field B (x) is composed of an ordered field, Bo(x), and a turbulent (random), 
zero-mean component, Bt (x), such that 



B(x) = Bo(x) + Bt(x). (4) 

For simplicity, we assume the polarized flux F (x) to be composed of an ordered component only (the 
more general case where a turbulent component is added was discussed in Paper II). Stationarity, 
homogeneity and isotropy (see Sec. 15. 2p in the magnetic field strength were assumed for Equation 
([U, while statistical independence between ordered and turbulent components will also be implied 
in what follows. 

Following the analysis of Paper II it is found that the dispersion function 1 — (cos [A^ (i)]) can 
be broken up into turbulent and ordered terms 



1- (cos[A$(^)]) = [b'^ (0) - b'^ {£)] + [a'^ (0) - a'^ {£)] 

= {b^O) + [a^O)-a^i)]}-b^£), (5) 

with the ordered and turbulent autocorrelation functions given by 
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^ ^ (B-B(0)> ^ ^ 

^ ^ (B.B(0)> ' 

respectively. The ordered function (0) — o? {£) can be advantageously modeled with a Taylor 
series 



oo 

a^{0)-a^{e) = Y,a2ji'^^, (8) 

while 6^ (0) is simply the turbulent to total magnetic energy ratio (that is, for the corresponding 
cloud- and beam-integrated quantities). As was shown in Paper I and Paper II the quantity within 
curly braces in Equation ([5]) can be readily determined from polarization maps by calculating the 
dispersion function (i.e., the left-hand side of that equation) from the data and fitting 



oo 

6^ (0) + [a^ (0) - {£)] = (0) + as,^'^' (9) 

to the data outside of the region where 6^ {i) is dominant (i.e., at lower values of t). Once Equation 
([9]) is evaluated the (normalized) turbulent cloud- and beam-integrated autocorrelation function 
h"^ {I) can be extracted from the data through Equation ([5]). 

Alternatively, the integrated turbulent autocorrelation function 6^ (£) can also be analytically 
derived using (see Eq. [A5] of Paper II) 



(BfBt(£)) 



H (a) H (a' + t) 



2 I " / u , 

1 - — I 7^3D,t [v, u) du 



A 



A 



d^a'd^a, (10) 



with 7^3D,t {v, u) = {F (a, z) F (a', z')) (Bt (a, z) ■ Bt (a', z')), u = \z' — z\ and v = |a' — a|, and a 
similar equation for {^B-B(O)). Since we are mostly interested in determining the shape of the 
magnetized turbulence power spectrum, we will concentrate on the Fourier transform of iP' {£) (see 
Eq. [A 12] of Paper II) 



1 



\H{K 



dku 



(11) 
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where ^-^^^ = (B-B (0)) and the Fourier transform of a function is represented by simply replacing 
the spatial arguments by their k-space counterparts (e.g., 7^3D,t iv,u) ^ 7^3D,t O^v, ku))- Equation 
pip can be redefined with 

b^K) = \\H{K)f'^^, (12) 
B 



where TZt (k^,) = J 7^3D,t O^v, ku) sinc^ (fc„A/2) dk^ is now interpreted as the two-dimensional turbu- 
lence power spectrum we seek to evaluate. We will accomplish this by taking the Fourier transform 
of 6^ {£) obtained from the data, as explained in the discussion following Equation ([9]) above, to 
evaluate (k„) on the left-hand side of Equation (jl2p and then invert this relation to determine 



7^t (k^,) / (b"^) (with a Wiener filter to remove the filtering due to \\H (kt,)|| ; see Appendix). 



Observations 



The observ ations for Orion KL were carried out on 10 September 2006 and 6 January 2008 
using the SMA ()Ho et al.l[20o3 lN in the compact array configuration, with the projected baseline 
lengths ranging from 15 to 80 kA (A = 870 /im). The phase center is at RA(J2000) = 5^35™14'!5, 
Decl(J2000) = — 5°22'30'.'4. The SMA receiver system has two sidebands, each with a bandwidth of 
~ 2 GHz. The sampled sky frequencies range from 345.5 to 347.5 GHz in the upper sideband and 
from 335.5 to 337.5 GHz in the lower, with a uniform spectral resolution of 0.812 MHz (corresponding 
to a velocity resolution of 0.7 km s~^). At these frequencies, the primary beam size (or field of view) 
of the SMA is ~ 32". Within the observational bandwidth, there is a significant contribution to the 
total emission from spectral lines of a number of molecular transitions (notably CO (7 = 3—7-2) and 
SiO (J = 8 — ?• 7)), and the continuum is generated after removing the spectral line contamination. 
Using natural weighting of the visibilities, the synthesized beam size is 2'.'6 x 1'.'7. The noise level 
in the Stokes I image is ~ 0.3 Jy beam~^. This is much higher than the theoretical noise level due 
to the limited dynamic range in the Stokes I map. The noise levels of both the Stokes Q and U 
images, which are much fainter, are much closer to the theoretical noise level, at 10 mJy beam~^. 
Our observations for Orion KL have much da. ta in common with those used for the polarization map 
presented in Figure 2 c) of iTang et al.l (|2010l ). which we refer the reader to in view of its similarity 
to the map that can be derived from our data. We note, however , that our data h av e a slightly 
higher spatial resolution than the 2'.'8 x 1'.'8 synthesi zed beam si z e of iTang et al.l (j2010l ). iKoch et al. 



(j2010l ) also performed a dispersion analysis on the iTang et al.l (j2010l ) map based on our Papers I 
and II. Importantly for our analysis, we use all polarization vectors available at the sampling rate 
of 0'.'25, provided they satisfy the condition p > 3ap, with p and dp the polarization level and its 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 
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uncertainty, respectively. This implies that many of our data points are correlated with each other 
since several are contained with a single synthesized beam profile. This is beneficial for our analysis 
as it allows for a better fit of Equation ([9]) from the dispersion data. The correlation between data 
points is also accounted for in the evaluation of the dispersion function and its uncertainty as a 
function of I (see Appendix B of Paper II). 



The data were calibrated and processed using the software package MIRIAD (jWright fc Sault 



19931 ). The gain calibration was obtained from observations of the QSO 0528+134. It is necessary 
to remove the contributions due to instrumental polarization as these are roughly similar in m ag- 
nitude to the observe d source polarization and can corrupt the data (see iMarrone et al.ll2006l and 
Marrone &: Radl2008l for the details of this method). The instrumental polarization was obtained 
from observations of the strong quasar 3c273 for 2 hours during transit. The total intensity (Stokes 
/) map was deconvolved using the task CLEAN in MIRIAD. We derived the polarized intensity 
(Jp) and position angles (P.A.s) with the task IMPOL, also using the CLEANed Stokes Q and U 
maps. The task IMPOL further removed the effects of the bias of the positive measure of /p. 

The data for IRAS 16293 and NGC 1333 IRAS 4A were also obtained at the SMA and reduced 



in a similar maii r ier. T hey were previously published and described in detail in iRao et al.l (j2009l ) 
and Girart et al. (j2006l ) . respectively. Their synthesized beam sizes and sampling rates are 3'.'1 x 2'.'0 
and 0'.'25 for IRAS 16293, and 1'.'6 x I'.'O and 0'.'2 for NGC 1333 IRAS 4A. 



4. Results 

4.1. Dispersion Functions from Interferometry Data - Orion KL 

Although the discussion that follows applies equally well to all th e data presented in this paper 



we will refer only to the Orion KL data (see map of Figure 2 c) of iTang et al.l 120101 ) in this sub- 
section in order to discuss the application of the dispersion analysis to interferometry data. One of 
our main goals here is to emphasize that there are important differences with results obtained with 
interferometry in relation to what one would get with single-dish data. 

We first note that because of the unavoidable 180-degree ambiguity when determining the 
orientation of the magnetic field from polarization data, one has to be careful when proceeding 
with a dispersion analysis of polarization angles. Indeed, any such analysis can only be successfully 
carried out in regions of polarization maps for which changes in polarization angles with position 
are sufficiently smooth to ensure, at least to a reasonable degree, that there are no possible reversals 
in the field direction. For this reason, we have excluded from our analysis a small isolated "clump" 
located at an offset of Aa ~ 8" and A(5 ~ 11" to the north of position of pea k intensity of Orio n 



KL (i.e., at a (J2000) ~ 5'"35°'14!9 and 5 ( J2000) ~ -5°22'17" in Figure 2 c) of lTang et al.ll2ninl ;i. 



Following the analysis done in Paper II, the dispersion function 1 — (cos [A<I> (£)]) and turbulent 
autocorrelation function 6^ (£) determined from the Orion KL polarization map are shown as a 
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function of the distance I in Figure [T] (symbols in top and bottom graphs, respectively). It is 
important to note that although these functions are only plotted for ^ > 0, the data is actually 
two-dimensional in nature and exhibits cylindrical symmetry (in the plane containing I about the 
£ = axis) because of the assumed isotropy of the dispersion function. With the prescription given 
in Paper I (or Paper II) one then fits the (sum of the) turbulent to total magnetic energy ratio \? (0) 
and the ordered magnetic field component X]j'^2j^^"'\ for which we use the lowest order polynomial 
that fits the data, with the broken curve shown in the top graph (see Equation [9]); the subtraction 
of the data from that curve would then yield }? {£) from which the analysis could proceed (see 
Equation [S]). Although this is perfectly adequate for single-dish data, an important fact needs to 
be emphasized when dealing with interferometry data. 

The aforementioned subtraction of the dispersion data from the broken curve in the top graph 
of Figure [T] leads to turbulent autocorrelation data that satisfy {£) > for most if not all values 
of i, which also implies, of course, that the integral of that function over space will be positive. 
But this should not be possible with a polarization map obtained with an interferometer, as can be 
assessed from Equation ()12p . since 



27T / {£) £di = {K = 0) 
Jo 

1 







Hik,)\\'ntiK] 



k„=0 



(13) 
(14) 



because H (k^, = 0) = for the so-called "dirty beam" of an interferometer. This, of course, is 
directly related to the well-known missing-flux "issue" that is implicit to interferometry data. But 
importantly, it is also contrary to single-dish data where (k^, = 0) = 1, when inefficiencies are 
accounted for. It therefore follows that an acceptable fit 6^ (0) + ^ja2j^^-' to our dispersion function 
interferometry data should satisfy the condition 



poo 

2tt / b'^{e)id£ = 0. (15) 
Jo 

Although we have restricted the condition of Equation (fT5|) to the turbulent component 6^ (^) ^ 
it also applies equally well to the ordered component and the total normalized autocorrelation 
function of the magnetic field given by Equation ([1]). The important point we need to acknowledge 
is that it is impossible to know exactly what correct fit for the turbulent to total magnetic energy 
ratio and ordered component (i.e., fe^ (q) _|_ Y^ja2ji'^-') applies to our Orion KL data (or any other 
data set) that will verify Equation (fT5|) . as there is an infinite number of combinations for b^ (0) 
and X^ja2j^^"' that will satisfy this condition and no way to discriminate between them. We must 
therefore accept this as a fundamental limitation to the analysis when dealing with interferometry 
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data alone. This problem can only be avoided if the polarization map includes data that fill the 
low- frequency portion of the spectrum (including = 0), i.e., when single-dish data are available, 
since the condition given in Equation (|15p will then not apply anymore and the degeneracy on the 
aforementioned fit will be lifted. 

On the other hand, it is perhaps reasonable to expect that our fit for (0) + X]j^2j^^"' (as shown 

in Fig. [1] for example) will affect the shape of the turbulence power spectrum TZt (k^,) / (^b"^"^, which 
we seek to determine from the dispersion function, mainly at low spatial frequencies. More precisely, 
the spectral content associated with the fit 6^ (0) + X^ja2j^^'' is to a large extent concentrated at low 
frequencies and will have a diminishing effect on our determination of the shape of the magnetized 
turbulence power spectrum when moving to higher frequencies. For the purpose of the present 
analysis, we will proceed by neglecting the condition stated in Equation (jlSp and perform our 
study as we normally would for single-dish data, and then discard the low frequency portion of 
the turbulent spectrum from our analysis (see below). That is, we adopt the fit shown in the top 
graph of Figure [1] (i.e., the broken curve) for Orion KL for our analysis. Accordingly, the associated 
turbulent autocorrelation function 6^ {£) is that shown in the bottom graph of Figure [1] (symbols). 
Also plotted is the radial profile of the "mean autocorrelated beam" (broken curve). This profile is 
obtained by first computing the autocorrelation of the synthesized beam, since it is this beam that 
is used for the polarization map and appears in the expression for the dispersion function (see Eqs. 
[T] and |10|). and then averaged azimuthally in the same manner as are the dispersion data. This 
represents the contribution of the synthesized beam to the (width of) the turbulent autocorrelation 
function 6^ {£). That is, this is what {£) would look like in the limit where the intrinsic turbulent 
correlation length were zero. A comparison of this autocorrelated beam profile to that of 6^ (£) 
clearly shows the significant contribution of the magnetized turbulence TZt (^) to the overall width 
and shape of 6^ {£). 

The top graph of Fig ure [2I shows the spectra associated with 6^ (-^) (i.e., (/c) with k — [k^^lj 
symbols) and the mean autocorrelated synthesized beam (i.e., \\H (/c)||^; broken curve) calculated by 
taking the Fourier transform of the corresponding functions shown in the bottom graph of Figure 
[T] (see Appendix). Also shown is the spectrum of the visibility data (or dirty beam; dot-broken 
curve), which is normalized to its peak level (using the scale on the right). We can now readily 
verify that the telescope dirty beam is responsible for the condition stated through Equation (|15p 
since it vanishes at A; = 0, as expected. The fact that our data for 6^ (k) are non-zero at and close to 
A; = is a reflection of the fact that, as we discussed above, we cannot precisely evaluate it from the 
dispersion function. Because of this we exclude the first two points of the spectrum at fc/27r = and 
0.5 arcsec"^ from our analysis and concentrate on the rest of the spectrum (i.e., higher frequencies). 

The bottom graph of Figure [2] presents the results of our analysis for Orion KL. First, Equation 
P2p is inverted using a simple Wiener optimal filter (see Appendix), which is possible since the beam 
profile ||//(A;)||^ is well characterized, to yield the magnetized turbulence power spectrum profile 
TZt (k) I \B) (again with k = |k^|; symbols). 
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The tur bulence power spectrum is not usually expressed with TZt {k), however, but rather with 
(lFrischlll995h 



7^K {k') = 47^fc'27^3D,t (A;') , 



(16) 



which is a "one-dimensional" representation of (the three-dimensional) 7^3D,t (k^j ^m) with k''^ = 
|kt, p + k'^. This "Kolmogorov-like" spectrum is particularly well-suited for cases of isotropic tur- 
bulence. It is unfortunately not possible with our data to recover T^k (k') from TZt (k) because of 
the integration over k^ in Equation (jllD that is inherent to the measurement process. On the other 
hand, it is still possible to define another "one-dimensional" power spectrum TZib (k) with 



7^lD (k) = 2TTknt (k) , 



(17) 



which is obtained from our data in a straightforward manner. We show the corresponding result 
for Orion KL with the solid curve in the bottom graph of Figure [2j Since it is also customary to 
parametrize the turbulence power spectrum in the inertial range with a power law, we also plot 
such a fit to T^iD (k) / show that the inertial range approximately scales with /j-(2-9±0.9)_ 

Although this result is consistent with theoretical expectations, the small number of spectral points 
available for the fit and our lack of knowledge concerning the precise shape of the ordered component 
to be subtracted from the dispersion function reduces the robustness of our determination for the 
spectral index (as is exemplified by the significant uncertainty on its value). 

Another important parameter that characterizes the magnetized turbulence power spectrum 
is the cutoff frequen cy fcAD. or length scale Aad = 27r/fcAr), at high frequeiicies, which is likely due to 
ambipolar diffusion ( Li &: HoudelboOS : Hezareh et al .112010 : Lazarian et al. I booi l iFalceta-Gongalves et al. 
2010l : lTillev fc Balsaral[201ol : lHoude et al.lboill ). We find from our results for TZm {k) / <^^^|) in Fig- 
ure [2] that k\D/2TT 0.2 arcse c~^. Although this value (Xab ^ 11 mpc, see below) is also consistent 



with theoretical expectations ([Lazarian et al.ll2004l : iFa 



2010 ) as well as observationally determined values ( Li fc Houdel 2008^ 



ceta-Goncalves et al 



201C: 



Tilley &: Balsara 



Hezareh et al 



201 



37 



we 



should ensure that this spectral cut off is real and not artificially imposed by data processing or 
the finite spatial resolution of the interferometer. This can be verified with Figure [3] where we have 
magnified the vertical scale of the top graph of Figure [2l It is observed that this spectral cut-off is 
present in the data for (k) (i.e., before applying the Wiener filter) and seen to happen well within 
the bandwidth subtended by the synthesized beam (or rather that of \\H (A:)||^), which cuts off at 
k/2TT ~ 0.6 arcsec"^. This plot also shows that kxD/2-K ~ 0.22 arcsec"^, or alternatively Aad — 9.9 
mpc for Orion KL (assumed to be at a distance of 450 pc). 



-10- 



4.2. IRAS 16293 



We applied our dispersion analysis to the polarization map of IRAS 16293 of lRao et al.l (|2009l ). 
The results for the dispersion and turbulent autocorrelation functions, the turbulent spectrum, 
and the determination of Aad are shown in Figures 31 [5l and [6l respectively. The dispersion and 
turbulent autocorrelation functions (Fig. [4]) display similar characteristics as those for Orion KL, 
but the differences are made more obvious when considering the turbulence power spectrum. Indeed, 
Figure [5] shows that the Kolmogorov-like spectrum for IRAS 16293 scales as ~ 

^-(1.8±0.3)^ while an 

inspection of Figure [6] makes it clear that the apparent cut-off in the spectrum at A; ~ 0.4 arcsec""^ 
(vertical dotted line) is likely due to beam filtering. We therefore only report an upper limit of 
•^AD ^ 1-8 mpc for IRAS 16293 (assumed to be at a distance of 150 pc) for the high frequency 
spectral cut-off due to ambipolar diffusion. 



4.3. NGC 1333 IRAS 4A 

The result s of ou r dispersion analysis as applied to the NGC 1333 IRAS 4A polarization map 



of lGirart et al.l ()2006l ) are shown in Figures [71 [U and[9l For this source the Kolmogorov-like power 
spectrum is observed to scale with ~ /j-(i-4±0.4) ^ggg pjg_ [gjj^ while we find kpj)/2T: ~ 0.66 arcsec"^ 
(or Aad — 2.2 mpc with an assumed distance of 300 pc) for the high frequency spectral cut- 
off due to turbulent ambipolar diffusion (vertical dotted line in Fig. [9]). However, because of the 
weakness of 6^ {k) about A;aD) the latter's proximity to the spectral cut-off due to beam filtering, and 
the aforementioned uncertainty in fitting the ordered component of the turbulent autocorrelation 
function, we must acknowledge that this estimate for Aad is tentative. As was the case for Orion 
KL (and IRAS 16293 for the spectral index) these parameters are consistent with expectations. 



5. Discussion 

5.1. Limitations of the Dispersion Technique 

The high spatial resolution with which the polarization data analyzed in this paper were ob- 
tained has allowed us to determine fundamental parameters that characterize the magnetized tur- 
bulence power spectrum in some well-known star-forming regions. However, the same observing 
mode that allows these realizations also brings with it some limitations due to the filtering of low 
spatial frequencies inherent to interferometry. We already discussed in detail in Section H] how this 
impedes the precise determination of the power spectrum at low frequencies. We now discuss two 
more consequences that result from this limitation. 
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5.1.1. The Chandrasekhar- Fermi Technique 

As was shown in Paper I, the value for the turbulent to ordered magnetic energy ratio {i = 0) 
obtained when fitting the dispersion function (see the broken curves in the top graphs of Figs. [H HI 
and [7]) corresponds to the value that one would normally use for t he determination of the ordere d 



magnetic field strength with the Chandrasekhar- Fermi equation (IChandrasekhar fc Fermil Il953l ) . 
More precisely, we have (see Eqs. [7] and [8] of Paper I) 

where p and a (v) are the mass density and the one-dimensional turbulent velocity dispersion, 
respectively. However, as was noted earlier the value of 6^ {i = 0) determined from interferometry 
data alone cannot be precisely determined because of the filtering of low spatial frequencies. An 
error in the estimate of the ordered magnetic field strength will then follow due to the presence 
of b'^ {i = 0) in the denominator of Equation (fT8|) . Although such an error will also arise with 
single-dish data (where the filtering happens instead at high frequencies), the relative importance 
of this error in comparison to what is expected with interferometry observations can be studied by 
considering the Fourier transform that links 6^ {£) to 6^ (k^,) 

6' (i) = II b' (K) e^-'^^d'K. (19) 

It follows from this and Equation (fT2|) that 



b" (^ = 0) = II b^ iK) 

= ^-j-^ ff \\H{K)\?'R,{\^,)<fK, (20) 

(27r)2(S^\ J J 



which is valid in general. For a given magnetized turbulence power spectrum TZt (k^,) the difference 
in 6^ (^ = 0) obtained with interferometry and single-dish observations resides in the nature of the 
filtering H (k„) applied to the data. Although a determination of b'^ {£ = 0) with single-dish will 
also be imprecise because of the spectral filtering at higher frequencies, the error is potentially more 
significant with interferometry since the spectral filtering is concentrated at low frequencies (see the 
"visibility" curves for Figs. O [5l and [8]) where the turbulent spectrum TZt (k^) peaks. Evidently, 
the relative importance in these errors will depend on the precise shape of the corresponding single- 
dish and interferometer dirty beams, as well as that of the underlying spectrum. Moreover, the 
nature of the imprecision in the evaluation of 6^ {£ = 0) is made more complicated by the fact that 
(b^/ = (B-B (0)), present in the denominator of Equation (|20|) . will also contain the same filtering 
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integral as well as a similar one for the ordered component of the magnetic field. In some cases this 
may alleviate the aforementioned error, in others it may worsen it. 

In the analysis of the SHARP OMC-1 data presented in Paper II the effects of the beam 
filtering and signal integration through the thickness of the cloud were corrected for by modeling 
the dispersion function while assuming circular Gaussian magnetized turbulent autocorrelation and 
beam functions. It is important to note, however, that even if it were possible to completely remove 
the filtering due to |li?(kt,)||^ in Equation (|20p . the corresponding value obtained for (i = 0) 
would still be underestimated because of the aforementioned signal integration along the line of 
sight. A correction for this effect would require a determination of the turbulent correlation length 
(see Paper II), which can be obtained by measuring the spectral width of TZt (k„)- More precisely, 
the turbulent correlation length is inversely proportional to the width of TZt O^v)- As this could only 
be achieved in general if TZt 0<-v) is known at low frequencies (where it peaks), it is apparent that 
one would greatly benefit from combining single-dish and interferometry observations to maximize 
the spectral coverage at both ends of the spectrum. 



An inspection of the magnetized turbulence power spectrum TZt {k) / \B J (or TZm {k) / {B j) 
shown in Figures [2l HI and [8] makes it clear that we are able to determine its shape only at the high 
frequency end, whereas it is expected that the inertial range of the power spectrum should extend 
over several decades in length scale (or spatial frequency) . Taking the case of Orion KL as an example 
(see Fig. [2]), and acknowledging the fact that, as discussed in Section [H our spectrum is unreliably 
estimated on the low frequency end (i.e., for k/2-K < 0.05 arcsec"^) we find that our analysis uncovers 
much less than a decade of the underlying spectrum (i.e., O.OTarcsec^^ < A;/27r < 0.2 arcsec"^). 
One must therefore be cautious in putting too much weight in our determination of the scaling laws 
characterizing the the small portion of the inertial range probed with our observations of Orion 
KL, IRAS 16293, and NGC 1333 IRAS 4A. Correspondingly, we once again emphasize the benefits 
that would thus be gained by combining single-dish and interferometry data. That is, a single-dish 
map of suitable spatial extent would ensure a good low frequency coverage, while a high-resolution 
interferometry map of the same region would extend the measured turbulence power spectrum far 
enough to precisely characterize the inertial range as well as the turbulent ambipolar diffusion scale. 

We should also keep in mind that the three-dimensional turbulence power spectrum underlying 
our data in Equation (|10p . i.e., i?3D,t {v-, u) = {F (a, z) F (a', z')) (Bt (a, z) Bt (a', z')), contains the 
autocorrelation of the (ordered) polarized emission as well as that of the turbulent magnetic field. 
It therefore follows that the magnetized turbulence power spectrum we extract from our data does 
not exactly correspond to that of the turbulent magnetic field, but is somewhat broadened by the 
polarized emission spectrum. The importance of this effect may be advantageously investigated 
through numerical analyses and simulations. 



5.1.2. The Magnetized Turbulence Power Spectrum 
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5.2. Further Improvements and Applications 



In all of our applications of the dispersion technique (i.e., in Paper I, II, and here) we always 
treated turbulence as being isotropic with a Kolmogorov-like power spectrum. This is, of course, a 
simplification that we do not expect to hold for magnetized turbulence in a weakly ionized plasma 
due to the anisotropy of motions in directions parallel and perpendicular t o the magnetic field 
broug ht about by the Lorentz force. If for example one considers the theory of lGoldreich &: Sridhar 
(jl995l ) for incompressible magnetized turbulence, then different power law scaling are expected 
for the power spectra meas ured along these tw o distinct orientations. Such anisotropy has, in fact, 
been measured in Taurus by lHeyer et al.l ( 20081 ) through ^^CO ( J = 1 ^ 0) observations and optical 
polarization measurements using principal component analysis. 

It would be straightforward in principle to extend our dispersion technique to allow for the 
detection of anisotropy. We would simply have to locally determine the mean orientation of the 
magnetic field about the position of a given datum on a polarization map and define, say, two sets of 
displacements i± and £|| depending whether the distance vector i linking that point to another one 
on the map is oriented more or less perpendicular or parallel to the mean magnetic field, respectively. 
Two dispersion analyses could then be performed, one for each of the i± and data sets. The main 
constraint in applying this technique for the polarization maps analyzed in this paper is the lack 
of data points. One needs a large number of points in order to accurately estimate the dispersion 
function. Moreover, it is also imperative that the change in orientation of the polarization vectors 
on a map varies smoothly enough that a mean direction for the magnetic field can be adequately 
calculated. Accordingly, we plan to attempt the implementation of this technique on the SHARP 
OMC-1 polarization map presented in Paper II in a future publication. 

Another natural extension of this technique concerns th e analysis of po l arizat ion maps of face- 
on spiral galaxies (see, e.g., the map of M51 presented in iFletcher et al.l (j20ld )). Although the 
polarization measured for external galaxies is not due to emission from anisotropic grains but from 
synchrotron radiation, we see no reason why the dispersion technique could not be applied to 
such cases. It would then also be natural to study the aforementioned anisotropy of magnetized 
turbulence since face-o n spiral galaxies oft en show a mean magnetic field orientation that closely 
traces the spiral arms ([Fletcher et al.ll2010l ). Such analyses would also provide a detailed study of 
magnetized turbulence on a much larger scale than we have achieved so far (i.e., on galactic scales 
instead of that of molecular clouds). 



6. Summary 

We presented an application of our magnetic field dispersion technique to high spatial resolution 
SMA polarization maps obtained for Orion KL in OMC-1, IRAS 16293, and NGC 1333 IRAS 4A. 
We showed how one can take advantage of such high resolution data to characterize the magnetized 
turbulence power spectrum in the inertial and dissipation ranges. For Orion KL we determine that 
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the inertial range of the spectrum approximately scales with /;;-(2-9±0.9) report a value of 9.9 

mpc for AaD) the high spatial frequency cutoff presumably due to turbulent ambipolar diffusion. 
For the same parameters we have ~ and a tentative value of Aad — 2.2 mpc for NGC 

1333 IRAS 4A, and ~ A:-(i-8±o.3) ^^^^ ^ppgj. ji^^^ ^^^^ < i g ^j. jj^^g ig293. 



The authors thank the referee, P. M. Koch, for his valuable comments, which greatly improved 
this paper. M.H.'s research is funded through the NSERC Discovery Grant, Canada Research Chair, 
Canada Foundation for Innovation, Ontario Innovation Trust, and Western's Academic Development 
Fund programs. 

A. Wiener Filter and Data Processing 

If we take into account the contribution of the spectral noise n (kt,) to 6^ (k^,) we can write 



iK)=P (k,)+n(k,). 



(Al) 



where the function 



(K) = \\H{k, 



B 



(A2) 



is assumed noiseless. The well-known solution for the Wiener filter to be applied to 5^ (k„) in order 
recover TZt (k^) / \ Bj (in the least-squared sense) is given by ( Press et al. 1992) 



■(k.) 



P (k.) 


' l\\H{K)\? 


b^ (k.) 


' + ||n(k,)f 



(A3) 



It is then necessary to somehow estimate the noise level n (k^) and insert Equation ()A2p into 
Equation ()A3P to express the filter as a function of ||i?(kt,)||^ and the signal-to-noise ratio, as is 
often done. It is straightforward to show that the Wiener filter for our problem can be expressed as 



\H{K) 



+ 



\\n{V.v) 



\n,{v.)f/{B'') 



(A4) 



Equation ()A4p is seen to tend to the obvious limit of (k^,) = 1/ \\H (k^,) || when n (k„) vanishes. We 
compute our Wiener filter by i) determining the mean level n for n (k^,) in the high frequency end of 
the spectrum for 6^ (k^,) where tp' (k^,) is negligible, ii) subtracting n from 6^ (k^) to approximately 
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2 \ 2 2 2 / 2 N 

obtain TZt{}i.v) / \B ), and in) inserting ||n|| (in lieu of ||n(k„)|| ) and ||7^t(ki,)|| / \B 
Equation ([X4]) . 



m 



Finally, we note that we processed the data in the spectral domain using discrete Fourier 
transforms (DFT). Although the shapes of the Ip' {€) functions shown in Figures [1] |4l and [7] are 
such that they tend to smoothly approach zero at the larger values of ^, the presence of noise and 
of residual levels can potentially be the cause of "edge effects" and ensuing contamination of the 
spectra 6^ (kt,). We have therefore windo wed the data and the synthesized beam in the £ domain 
with a Hanning window (jHammine] 119971 ) before applying DFTs in order to minimize these effects. 
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Fig. 1. — Top: The dispersion function 1 — (cos [A$ (i)]) for Orion KL. The broken curve ("ordered") 
is the fit for the sum of the turbulent to total magnetic energy ratio 6^ (0) and the ordered component 
Yl^=i o^2j^'^^ implicit to the data (symbols) using the points where i > 1" . Both functions are plotted 
as a function of i. Bottom: the turbulent autocorrelation function }? {£) (symbols), as obtained by 
subtracting the data points to the "ordered" curve in the top graph, while the broken curve shows 
the radial profile of the "mean autocorrelated synthesized beam" (see text). 
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Fig. 2. — Top: spectra associated to b'^ (i) (i.e., b'^ (k) with k = |kt,|; symbols) and the mean 
autocorrelated synthesized beam (i.e., ||i?(A:)||^) calculated by taking the Fourier transform of 
the corresponding functions shown in the bottom graph of Figure [TJ The visibility data (or the 
spectrum profile of the dirty beam) is also shown (see text) and is normalized to its peak level 
(using the scale on the right). Bottom: our results for the magnetized turbulence power spectrum 
profile TZt (k) / ^-^^^ (symbols) and the associated one-dimensional "Kolmogorov-like" spectrum 

7^lD (k) I {b^)- The broken curve shows an approximate power law fit k (2-9±0-9) to the inertial 
range. 
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Fig. 3. — Magnification of the vertical scale of the top graph of Figure [2j For Orion KL we detect 
kxT) /2t: ~ 0.22 arcsec^^ (or Aad — 9.9 mpc with an assumed distance of 450 pc) for the high 
frequency spectral cut-off, which is probably due to turbulent ambipolar diffusion (vertical dotted 
line). 
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Fig. 4. — Same as Figure [T] but for IRAS 16293; data points where £ > 6.25" were used for the fit 
on the top graph (broken curve). 




Fig. 5. — Same as Figure [2] but for IRAS 16293. However, for this source the broken curve on 
bottom graph shows an approximate power law fit /i;-(i-8±0-3) to the inertial range. 
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Fig. 6. — Same as Figure [3] but for IRAS 16293. However, for this source we observe that the high 
frequency cut off present in the data for 6^ [k) (i.e., before applying the Wiener filter) is correlated 
with the bandwidth subtended by the synthesized beam||ff (A;)||^, which cuts off at k/27r ~ 0.41 
arcsec"^. We therefore have an upper limit of Aad ^ 1-8 mpc for IRAS 16293 (assumed distance of 
150 pc) for the high frequency spectral cut-off expected from turbulent ambipolar diffusion (vertical 
dotted line). 
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Fig. 7. — Same as Figure [T] but for NGC 1333 IRAS 4A. The broken curve on the top graph 
("ordered") is the fit for the turbulent to total magnetic energy ratio (0) and the ordered component 
X^j=i o,2ji'^-' implicit to the data (symbols) using the points where £ > 3'.'4. 
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Fig. 8. — Same as Figure [2] but for NGC 1333 IRAS 4A. However, for this source the broken curve 
on the bottom graph shows an approximate power law fit A;-(i-4±0.4) ^.j^g inertial range. 
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Fig. 9. — Same as Figure [3] but for NGC 1333 IRAS 4A. However, for this source we tentatively 
detect A:AD/2vr ~ 0.66 arcsec"^ (or Aad — 2.2 mpc with an assumed distance of 300 pc) for the high 
frequency spectral cut-off, which is probably due to turbulent ambipolar diffusion (vertical dotted 
line). 



